Ring exchanges and the supersolid phase of 4 He 
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Using Path Integral Monte Carlo we calculate exchange frequencies in bulk hep 4 He as atoms 
undergo ring exchange. We fit the frequencies to a lattice model and examine whether such atoms 
could become a supersolid, that is have a non-classical rotational inertia. We find that the scaling 
with respect to the number of exchanging atoms is such that superfluid behavior will not be observed 
in a perfect 4 He crystal. 
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Recent torsional-oscillator observations by E. Kim and 
M. H. W. Chan on solid 4 He, both in the disordered 
absorbed vycor[l| and in bulk 4 He 0], have revived in- 
terest in the super-solid phase. In this phase, one has 
both long range crystalline order and superfluidity. Ex- 
periments find superfluid response to rotation at tem- 
peratures below 0.2K. The low temperature superfluid 
density is about 1% of the solid density and is found 
to be independent of density within a pressure range of 
25 to 60Bars. In the following, we assume the observed 
phenomena is described by equilibrium thermodynam- 
ics of pure bulk helium. In practice, whether such is 
the case, is one of the most important issues. Assuming 
equilibrium 0], a supersolid is characterized by both long- 
range translational order and by a non-classical response 
to rotation (NCRI). 

In proposing such a state, Andreev and Lifshitzjj] pos- 
tulated the existence of suitable defects (e.g. vacancies) 
and then estimated the properties of the dilute system 
of bosonic defects. The superfluid transition tempera- 
ture of point bosons is T c = 3.31p 2 / 3 7ir/m where p is 
the density of the bose condensing quasi-particles and 
m their effective mass. Taking the measured value for 
T c = 0.2K , writing the mass in terms of the bare helium 
atom mass m = fxmi and their density as a fraction c of 
the bulk helium density at melting p = 0.029A -3 we find 
that c ~ 0.012/^ 3 / 2 : if the quasi-particle is heavier than 
a helium atom, at least one out of every 83 lattice sites 
must have a defect; if the mass is 1/10 the atom mass, 
the relative density of defects will need to be 4 x 10~ 4 . 
Such a dense collection of defects could not be associated 
with the much lower concentration of 3 He impurities. 



Experiments, particularly the NMR experiments . 
on solid 3 He rule out zero point vacancy concentrations 
of more than 1 part in 10 14 . Pederiva et al.Q have done 
calculations of the energy of vacancies and find an energy 
of about 15 K±4K in the hep phase of bulk 4 He at the 
melting density, and 30 K at 50 bars, in agreement with 
various experiments, most notably X-Ray scattering^. 
Using path integral calculations (described below) we es- 
timate the energy to create an interstitial at the melting 



density is 48K± 5K; they are even more improbable than 
vacancies. Hence, the energies of point defects are sub- 
stantially greater than zero; in equilibrium there should 
not be enough of them at 0.2 K to bose condense; this 
precludes the Andreev mechanism for a supersolid. 

In this paper we examine the possibility that bulk solid 
hep helium, assumed to be free of defects such as impu- 
rities and vacancies, could have a supersolid phase. One 
might think that there would always be ground state de- 
fects, arising from the large quantum zero point fluctu- 
ations. Near melting, the r.m.s. vibration about the 
lattice site is 30%, so that at any instant of time, a good 
fraction of atoms are closer to a neighboring site than to 
their home site. However, the absence of an atom from 
a lattice site is not sufficient for having a supersolid; if 
the empty site is always accompanied by doubly occupied 
site, there can no mass current. Chester 9] proved that 
any Jastrow (i.e. pair product) wavefunction of finite 
range has both BEC and vacancies. However, Jastrow 
wavefunctions crystallize with much difficulty and with a 
transition density off by an order of magnitude 10] so that 
they are not reliable enough to be used to predict super- 
fluidity. In general, variational wavefunction approaches 
are suspect, since the energy is not sensitive to the low 
probability regions of configuration space that are impor- 
tant for the superfluid density and bose condensation. 

LeggettQ has written an upper bound to the super- 
fluid density in terms of the zero temperature density 
of the solid. Using densities from path integral calcula- 
tions, we find that the Leggett bound gives p s /p < 0.16 at 
the melting density, a value consistent with experiment. 
In this paper, we calculate values for ring exchange fre- 
quencies and using them find that the superfluid density 
should be zero in solid helium. 

To determine whether bulk helium could be supersolid, 
path integrals give a much cleaner frame-work; they can 
be used to compute the superfluid density and the mo- 
mentum distribution without the assumption of a trial 
wave function or any other uncontrolled approximation. 
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The partition function of N bosons is: 

Z = JnzZ I dR{R\e- 0H \PR) (1) 
• P J 

where H is the Hamiltonian, (3 the inverse temperature 
and R = {i"i, r 2 , . . . rjy}. For numerical calculations, the 
density matrix operator is expanded into a path, begin- 
ning at the configuration R and ending at PR. In terms 
of these paths, the superfluid density {i.e. the number of 
atoms not moving with the walls of the torsional oscilla- 
tor) is given by: 

^ = m < WW > 

p ~ h 2 f3N [ ' 

where W — L dtJ^iLi d,Yi(t)/dt is the winding number 
of the path around a torus. It is only exchanges on the 
order of the sample size that contribute to the superfluid 
density; local exchanges make no contribution. Using the 
PIMC method, superfluidity and freezing, happen natu- 
rally at the right density and temperature, without im- 
posing them in any way. The technical complications (ll) 
concern ergodicity of the random walk, and finite size 
effects: one has to take the limit as N oo. 

PIMC calculations find a superfluid density on the or- 
der of 3% at melting density (molar volume 21.04 cm 3 ) 
and about 1.2% at 55 bars (molar volume 19.01 cm 3 ) 
in a 48 atom (3x4x2) hep supercell. The superfluid 
density is larger than what is observed, and has a rela- 
tively weak pressure dependance. However, a cell with 
180 atoms, has zero superfluid density. We cannot be 
sure that the lack of winding paths has a physical origin, 
or is due to a lack of ergodicity within the PIMC random 
walk. To change from one winding configuration to an- 
other involves a global move of the paths which becomes 
very unlikely as the box size gets large. To avoid this 
problem, we turn to a PIMC approach which directly 
estimates individual exchange probabilities. 

The Thouless^^ theory of exchange in quantum crys- 
tals assumes that at low temperatures, the system will 
almost always be near one of the particular TV! arrange- 
ments of particles to lattice sites with rare, rapid, tun- 
ncllings from one arrangement to another. Hence, we 
label the particles with their initial lattice sites. Then 
the partition function is written as a sum over permuta- 
tions of sites onto themselves. We break up the permu- 
tation into cyclic exchanges \p\ 1 p2 ■ ■ -p n }', each cycle is 
independent, so that: 

n P 

z = z a J2HfpM- (3) 

P i=l 

where Zq = J dR{R\e^ l3H \R) is an uninteresting phonon 
partition function at low temperature. The contribution 
for a cycle: 

fvifi) = y J dR(R\e- fjH \ P R) = J P f3. (4) 



is proportional to (3 because it is localized in time; the 
coefficient, J p is the exchange frequency. 

A special case of the above par tition function is the 
Feynman-Kikuchi (FK) model[13j that assumes that 
Jp = Joe~ aLp where L p is the cycle length. The par- 
tition function is: 

Z = Z ^(/3Jo)' if 'e- Q S"" 1 L «. (5) 
p 

We make contact with previous work by assuming [3 Jo = 
1. Calculations were done assuming a cubic 3D lattice, 
first by Feynman and Kikuchi analytically and then nu- 
merically by Elser^4|- It i s found that the superfluid 
transition occurs as a becomes smaller from a localized 
state (small cycles) to large cycles for a c w 1.44. To un- 
derstand this critical value, consider the free energy of 
adding a link to an existing loop. A cubic lattice has 
coordination number 6 but a return to the previously 
visited site is impossible since it is a cyclic permutation. 
Adding a link costs probability e~ a but the entropy of 
the new link is 5; hence a c ~ ln(5) = 1.6. The critical 
value is smaller because of the 'self-avoiding' restriction 
within a cycle and on overlapping cycles. Using the same 
argument for the hep lattice gives a critical coupling of 
a c w 2.3. 

Using methods [rjlllfij developed for solid 3 He, we cal- 
culate the exchange frequencies and estimate how close 
they are to the critical value. We assume the helium in- 
teract with a semi-empirical [l^j pair potential. The fre- 
quencies for 2, 3 and 4 atom exchanges is very small[l3, 
e. g. J2 ~ 3fiK at melting density. However, small cyclic 
exchanges are quite different from the long exchanges 
needed to get a supersolid. We have performed exchange 
calculations of 50 different exchanges involving from 5 to 
10 atoms. All exchanges involve nearest neighbors, since 
calculations show that next-nearest exchanges are much 
less probable. We obtain accuracies on the order of 5% 
for the 5-particle exchanges and 10% for the 9 particle ex- 
changes. More than half of the exchanges involve winding 
around the cell boundaries, important because they are 
representative of the type of exchanges in a supersolid. 

Fig. shows the results of calculations of the fre- 
quency of the simplest winding exchanges: straight line 
exchanges in the basal plane. As assumed in the FK 
model, we find that the exchange frequencies decrease 
exponentially with the length of the exchange with an 
exponent of a = 2.64 near the melting volume 21.04 cm 3 , 
and a = 3.14 at the volume 19.01 cm 3 corresponding to 
P ~ 60 bars. 

To construct a more realistic model than the FK 
model, we need to take into account more details of the 
geometry of the exchange than just the number of ex- 
changing atoms. We assume that it is the internal ge- 
ometry of the exchange that matters 0; the detailed 
arrangement of the neighboring spectator atoms is much 
less important. In particular, we assume the action of a 
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FIG. 1: The exchange frequencies (J in K) versus exchange 
length L p for straight line exchanges in the basal plane that 
wind around the periodic cell. The inset shows the lattice sites 
in the basal plane for L p = 6; the red atoms show the atoms 
midway through the exchange. The circles with error bars 
are the hi(J) at two molar volumes; 21.04 cm 3 (black, open 
circles, solid line) and 19.01 cm 3 (red, solid circles, dashed 
line). Multiple points at the same value of L p are from cells 
with different number of atoms in the directions perpendicular 
to the exchange direction. The lines are least squares fits. 

given n-cycle (the log of the exchange frequency) is the 
sum over the internal vertices of the exchange: 

p 

J p = J exp[- a(flfc)]. (6) 
fe=i 

Here Ok are angles between successive displacements in 
the exchange cycle. In an hep lattice there are 7 possible 
angles between two nearest neighbor displacement vec- 
tors, but — only occurs in the pair exchange. We 
determine the parameters by fitting to the PIMC ex- 
change frequencies J p . Good fits are obtained; the model 
predicts the exchanges frequencies with an accuracy of 
about 20%. The resulting coefficients are shown in Fig. 
121 Typical errors on the coefficients are ±0.03. Note that 
there is a strong preference for exchanges that proceed in 
a straight line versus ones which have sharp angles; for 
them the incoming and outgoing particles are more likely 
to collide. The fitting coefficients are seen to depend lin- 
early on cos 4 (0/2) . We find no preference for exchange in 
the basal plane beyond the effect induced by the selection 
of angles. 

We also verified this model by comparing with results 
of the direct PIMC calculations (where permutations are 
generated dynamically) mentioned above and obtained 
rough agreement with the winding numbers and permu- 
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FIG. 2: (Color online) The model parameters vs. cos 4 (#/2) 
at the melting density (open circles). The angle is the inte- 
rior angle of the vertex as shown. The fits were done using 
50 exchanges in a 180 atom cell. The solid red circle corre- 
sponds to a pressure 60 bars. The two lines show a critical 
supersolid Hamiltonian with (solid line) and without (dashed 
line) angular dependance. 



tation cycle distributions. The 2- and 3- particle ex- 
changes are about twice the prediction of model. On 
the other hand, the exchanges in the 180 atom cell are 
somewhat smaller than Thouless theory. We expect such 
disagreements because Eq. Eldoes not contain the effect 
of spectator atoms on the exchange process, and small 
exchanges have not been included in the fit. 

The activation coefficient of the exchange is Jo = 7.2 ± 
0.5-ftT. This is the order of magnitude but smaller than 
the vacancy energy. If it is related to vacancies, then Jo 
should be of the same order of the vacancy energy at all 
densities, but on the contrary, we do not find a significant 
density dependance. 

We now analyze the model. Because the probability 
of retracing is small, we initially neglect the self inter- 
sections. This leads to a diffusion problem on a lattice, 
with a large probability of continuing in the same di- 
rection. Since a displacement only depends on the previ- 
ous displacement, it is an un-normalized Markov process. 
(Because the hep lattice has a basis, we have to label 
the displacements consistently on the A and B planes so 
that the transition probabilities are independent of the 
plane.) The probability of a given displacement vector II; 
will approach a steady state after many steps, [i refers 
to one of the lattice directions) and satisfy the eigenvalue 
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equation: 

Al1 .-- ( ? ) 

i 

By symmetry 11, can only depend on whether the direc- 
tion is in the basal plane, ILi, or out of the basal plane, 
n 2 . Define the partial sum D^ j — £\ e~ ai i where 
are either in the basal plane (1) or out of it (2), e.g. D u 
is the sum of the probability of both vectors being in the 
basal plane). The eigenvalue A, is the solution to the 
secular equation: 

(Dii-A)(D22-A)-D? 2 = 0. (8) 

Putting in the model parameters, we find A = 0.303 ± 
0.005. The probability of having an exchange of length 
p will equal a prefactor times A p , hence, since A < 1, 
our PIMC results imply that solid 4 He will have only 
localized exchange and thus cannot be a supersolid. 

Let us consider how the neglect of the self- intersections 
affects the critical value. Qualitatively, self-intersections 
must decrease the probability of long exchange cycles. 
To estimate the effect quantitatively, we perform random 
walks on the hep lattice and only count non-intersecting 
walks. If we artificially change the model parameters 
a by subtracting 1.2, the model becomes critical when 
A w 1.03 as opposed to A = 1. This line is shown in 
figure 2. Our computed value of A is much less than 
what is needed to allow for a supersolid. We expect that 
the presence of other exchange cycles will further increase 
the critical value of A. 

In summary, PIMC-computed exchange frequencies for 
hep solid 4 He show that only localized exchanges will 
be present and thus should not exhibit the property of 
nonclassical rotational inertia. Based on other theoretical 
and experimental findings, we think it unlikely that the 
observed phenomena are due to vacancies, interstitials or 
3 He impurities. Hence, one must look for an explanation 



of the experiments elsewhere, either to non-equilibrium 
effects or more complicated lattice defects. 
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